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Higher harmonics of anisotropic flow (vn with n > 4) in heavy-ion collisions can be measured either 
with respect to their own plane, or with respect to a plane constructed using lower-order harmonics. 
We explain how such measurements are related to event-plane correlations. We show that CMS 
data on va and ve are compatible with ATLAS data on event-plane correlations. If one assumes 
that higher harmonics are the superposition of non-linear and linear responses, then the linear and 
non-linear parts can be isolated under fairly general assumptions. By combining analyses of higher 
harmonics with analyses of V 2 and V 3 , one can eliminate the uncertainty from initial conditions and 
define quantities that only involve nonlinear hydrodynamic response coefficients. Experimental data 
on V 4 , Vs and ve are in good agreement with hydrodynamic calculations. We argue that V 7 can be 
measured with respect to elliptic and triangular flow. We present predictions for V 7 versus centrality 
in Pb-Pb collisions at the LHC. 

PACS numbers: 25.75.Ld, 24.10.Nz 


I. INTRODUCTION 


In the last year or so, LHC and RHIC experiments have 
probed anisotropic flow [T] and its fluctuations HIS] to 
an unprecedented degree of precision [IHH]- These new 
analyses include in particular detailed analyses of higher 
Fourier harmonics (u4, U5, vq) and their correlations with 
lower harmonics {v2, V3). The scope of this paper is 
twofold. The first goal is to point out specihc relations 
between seemingly different observables found in the re¬ 
cent experimental literature, and to propose new observ¬ 
ables. The second goal is to show that measurements of 
higher harmonics can be combined with measurements 
of lower harmonics in a way that facilitates comparison 
with theory. As an illustration, recent experimental re¬ 
sults are compared with hydrodynamic calculations. 

The CMS Collaboration has measured V4 and ve with 
respect to their own direction, and with respect to the 
direction of elliptic flow V2 [1] (see also 0); on the other 
hand, the ATLAS Collaboration has measured a large 
number of event-plane correlations In Sec. [TTj we 
clarify the relation between these observables and show 
how they are related to one another. In particular, we 
show that CMS and ATLAS data on V4 and ve are com¬ 
patible. We explain how odd harmonics, such as U5 or 
V7, can also be analyzed with respect to the direction of 
lower harmonics. 

While recent experimental data have been compared 
to several theoretical models, either event-by-event hy¬ 
drodynamic calculations [9HT2] or tranport models |13j . 
these comparisons offer little insight into the physics of 
higher-order harmonics. In particular, theoretical calcu¬ 
lations depend strongly on the model of the initial den¬ 
sity profile, which has long been recognized as the main 
source of uncertainty in modeling anisotropic flow Ha- 
On the other hand, there are hints that the physics of 
higher-order harmonics should be simple: for instance, 
the ratio V4/(v2)^ [IH [16] is equal to | at high trans¬ 
verse momentum pT in ideal hydrodynamics. 


In hydrodynamics, higher-order harmonics are super¬ 
positions of linear and non-linear response terms niHin]. 
This is recalled in Sec. Ill We explain how the linear 
and nonlinear terms can be isolated under fairly gen¬ 
eral assumptions. We show how analyses of higher-order 
harmonics can be combined with analyses of lower-order 
harmonics (v2 and V3) to form quantities which do not 
involve the initial state. These quantities are compared 
with hydrodynamic calculations. 

In Sec. |IV[ we list a few predictions for higher-order 
harmonics; in particular, we predict the value of vy, mea¬ 
sured with respect to V2 and V3, as a function of central¬ 
ity. 


II. OBSERVABLES FOR HIGHER HARMONICS 

Anisotropic flow is an azimuthal (p) asymmetry of the 
single-particle distribution m 

= ^ E w 

n— — oo 

where Vn = Vn exp(zn'I'„) is the (complex) anisotropic 
flow coefficient in the nth harmonic, and = V*. Both 
the magnitude [5^ and phase [1113] of Vn fluctuate event 
to event. 

The simplest observable involving Vn is a plain rms 
average [HIIS]: 

= vEW, (2) 

where angular brackets denote an average over events. 
The notation Vn{'^n} has been used earlier to denote the 
value analyzed with the event-plane method |4]. How¬ 
ever, the event-plane method does does not quite measure 
the rms average [13]. Therefore it should be replaced by 
the scalar-product method m, which is recalled in Ap¬ 
pendix]^ Note that our Unjlkn} is the same quantity as 
Vn{‘2} in the notation of the cumulant analysis |28j . 
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Alternatively, V4 can be analyzed with respect to the 
direction of V 2 [SSI ED], and Ve can be analyzed with 
respect to the direction of V 2 or that of V 3 , 


W4{^'2} 


Re(V4(V2*y} 

^/W^) 

R<v^{vi?) 

y/W^) 


(3) 


The triangular inequality implies |f4{5'2}| < U4{'1'4}, 
k 6 {^ 2 }| < 'Celt's}, kelp's}! < ^'6{^6}, i-e., V 4 and vq 
are larger when measured with respect to their own plane 
than with respect to another plane. The ratio of Vn{'^m} 
and (where n is a multiple of m) can be writ¬ 

ten as the Pearson correlation coefficient between Vn and 
which we denote by Pmn- 


^ Reiysq)^ _V4{^ 

Re(P6(^?)") t^6{^2} 

Re(P6C^3*)^) ^6{^3} 

\/(|V6P)(|R3|4) 


The correlations between event planes measured by 
ATLAS, which are denoted by (cos( 4($2 — 
(cos( 6($2 - 4'6))).u, and (cos( 6($3 - d> 6 ))).^, in Ref. [5|, 
are precisely P 24 , P 26 and p^e [13 ED- Note that the 
terminology “event-plane correlations” applied to such 
measurements is somewhat misleading, in the sense that 
these observables involve not only the angles of Vn, but 
also their magnitudes Ell- 

Figure [^presents a test of the first two lines of Eq. Q, 
where the left-hand side uses ATLAS data and the right- 
hand side CMS data. The overall agreement is very 
good, which shows that CMS and ATLAS data are com¬ 
patible, even though they are measured with different 
cuts of tranverse momentum px- Note that CMS uses 
the event-plane method, instead of the scalar-product 
method. This method yields a slightly lower correlation 
when the resolution is large [5|- This explains, at least 
qualitatively, why CMS data are slightly lower than AT¬ 
LAS data for midcentral collisions in Fig. 0 (a). 

While Pearson correlation coefficients are typically an¬ 
alyzed by integrating over all particles in a reference 
detector [ 3 , analyses of Vn with respect to a specific 
direction (either 4 '2 or can be done differentially, 
as a function of transverse momentum px [30] (see Ap- 
pendixj^for analysis details). Hydrodynamics predicts a 
slightly different px dependence depending on the refer¬ 
ence direction m- It is therefore interesting to general¬ 
ize Eq. (|^ to odd harmonics. V 5 and E 7 can be analyzed 
with respect to the directions of V 2 and V 3 in the follow¬ 
ing way: 


1'5{'I'23} 


Re{V^V^Vi) 

vwwm 
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FIG. 1. (Color online) Test of Eqs. Shaded bands corre¬ 
spond to the left-hand side measured by ATLAS [3 in Pb-Pb 
collisions at 2.76 TeV. Full circles correspond to the right- 
hand side, obtained using CMS data (Jj. 


^^7{df23} 


Re{Vj{Vi)^Vi) 

^/{\V 2 mV) 


(5) 


Quantitative predictions for these quantities will be pre¬ 
sented in Sec. |IV[ These projected harmonics are smaller 
than those defined by Eq.(0), namely, |r’5{4'23}| < 
UsilE's}, (and |?^7{4'23}| < r’rj'I'y}). The ratio of 
|w 5 {'I' 23 }| and usjTs} is again the Pearson correlation 
coefficient between V 5 and V 2 V 3 : 


Re(P5^2*^3*) _ ^ 5 (^ 23 } 


( 6 ) 


This quantity is very similar to the corresponding three- 
plane correlation measured by ATLAS [Sj: 


(cos(2$2 -I- 3$3 - 54>5))u, 


Re(P5^2*^3*) 

vrnFmmW) 


■ (7) 


More precisely, they coincide if the magnitudes of 
V 2 and V 3 are uncorrelated 0 namely, (|V 2 p|Vi 3 p) = 


^ A slight anticorrelation between IF 2 P and IV 3 P has been pre- 
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(|^2p)(|V3p). Throughout this paper, we use (cos(2$2 + 
3$3 — 5$5 ))uj from ATLAS as an approximation for p235- 
Note that even though U4{'I'2} and U6{4'2} are smaller 
than U4{5'4} and uel'I'e}; respectively, they are measured 
with better relative precision [J]. The reason is that these 
measurements use elliptic flow as a reference, which is 
measured very accurately. Triangular flow, ua, is also 
precisely known. We therefore expect that U5{'I'23} be 
determined with better relative accuracy than i'5{'l'5}. 
In the same way, we expect that even though no experi¬ 
ment has yet been able to detect a nonzero LHC 

experiments could already measure f7{'I'23}- 


Glauber simulations [ 18 ) . The idea behind this assump¬ 
tion is that V4L is produced by initial fluctuations in the 
fourth harmonic, which are not correlated with the mean 
eccentricity. Then, the last term in the right-hand side 
of Eq. Q vanishes, and the rms value of the linear part 
is 

V4L = V (|V4Lp) = \/ (u 4{'I'4})2 — (w4{'I'2})^ 

V5L = ^/W^) = - (^^5{'I'23})2- (10) 

This quantity has been measured as a function of cen¬ 
trality by the ATLAS collaboration [ 40 ) . 


B. Nonlinear response 

III. LINEAR AND NONLINEAR RESPONSE 


In hydrodynamics, anisotropic flow is the response to 
anisotropy in the initial density profile [ 34 ) . Harmonics 
V4 and higher can arise from initial anisotropies in the 
same harmonic ( 3 ) [ 3 HVf 37 ] (linear response) or can be in¬ 
duced by lower-order harmonics [HI | 38 l [ 39 ] (nonlinear 
response). To a good approximation | 5 D], one can write 

V 4 = V 4 L + X 4 {V 2 f 
V5 = V5L + X5V2V3 
14 = V(,L + X 62 (V 2 )^ -f 

1^7 = VtL + X 7 {^ 2 )'^V 3 , (8) 

where VnL denotes the part of due to linear response, 
and we have included the nonlinear terms involving the 
largest flow harmonics, V2 and V3. The interest of this 
decomposition is that the nonlinear response coefhcients 
X are independent of the initial density profile in a given 
centrality class [ 15 ]. We now explain how the linear and 
nonlinear parts can be isolated. 


A. Linear response 


The linear part of V4 and V5 can be isolated [ 40 ) by 
combining the observables introduced in Sec. |T^ Using 
Eqs. ([^ and one obtains 

(^4{'I^4})^ - (7^4{'I^2})^ = {\V4Ln - 
(^' 5 {'P 5 })" - (^ 5 {^ 23 })" = (I^ 5 LP) - '^^^^^^’' ( 9 ) 


These results are general: these combinations always sub¬ 
tract the nonlinear response. 

Erom now on, we further assume that the terms ap¬ 
pearing in the right-hand side of Eq. Q are uncorrelated. 
That is, we neglect the small correlation between the lin¬ 
ear and nonlinear parts which is seen in Monte-Carlo 


The nonlinear parts are obtained by projecting Eq. ([^ 
onto lower harmonics. Assuming again that the terms 
in the right-hand side of Eq. (|^ are uncorrelated, one 
obtains the following expressions for nonlinear response 
coefficients: 



W4{4'2} 

mi*) 

vmm 

{V3V*V*) 

1^5{'f'23} 

{\V2?\Vi\) 


{vdv^f) 

1^6{^2} 



{ VeiVif ) 

W{^3} 

mi*) 


( VrivmVi ) 

U7{4'23} 

(1^211^321) 



( 11 ) 


The left-hand side of these expressions can be calculated 
in hydrodynamics, and is independent of the model of ini¬ 
tial conditions, while the right-hand side can be inferred 
from experimental data. Eq. ( 11 ) therefore offers a direct 
comparison between hydrodynamics and data, where all 
dependence on initial state is eliminated jH]. Nonlin¬ 
ear response coefficients have been obtained using event- 
shape engineering [ 40 ] by the ATLAS collaboration. The 
present method does not require event-shape engineer¬ 
ing. The comparison between hydrodynamics and data 
is shown in Eig.[^ and we now explain in detail how these 
results are obtained. 

The numerators in the right-hand side of Eq. ( | 11 [ ) are 
the projected harmonics defined by Eqs. ^ and We 
use U4{'L2} and U6{'I'2} measured by CMS [ 1 ]|^ U5{'I'23} 
and are not measured directly, but can be in¬ 
ferred from P235 and pse using Eqs. ffl 

and Q. We use CMS data [ 1 ] for 1)5 and 
and ATLAS data [S] for P235 and p^Q. These correlation 
coefhcients, however, are expected to depend little on the 
experimental setup, as illustrated in Fig. [^ 


dieted in AMPT simulations |3m3,S| . but it is at most at the 10% ^ Since CMS uses the event-plane method, the results are slightly 

level. lower than the nominal quantities in Eqs. i3| mi- 
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FIG. 2. (Color online) Nonlinear response coefficients defined by Eq. (Ill as a function of centrality. Each panel corresponds to 
a different line of Eq. ( 111 . Dashed lines: ideal hydrodynamics. Solid lines: viscous hydrodynamics with 77 /s = 0.08. Symbols: 
experimental data (see text for details). 


The denominators in the right-hand side of Eq. ( 11 ) in¬ 
volve various even moments of the distribution of V2 and 
V3. There is no direct measurement of these moments 
to date. A straightforward procedure to analyze them 
is outlined in Ref. m- Alternatively, moments of the 
form can be inferred from cumulants [ 5 S]. The 

expressions of the hrst moments in terms of cumulants 
are: 


m^) = v2{2r 

{\V^f)=2v2{2}^-V2HV 

(|K|6) = 47;„{6}® - 9t-„{4}V{2}2 + 6u„{2}®. (12) 


For the moments involving both V2 and V3 (second and 
fourth line of Eq. (dll)), we further assume that the mag¬ 
nitudes of V2 and V3 are uncorrelated. 

Since different experiments have different acceptance 
(in particular in transverse momentum pt), it is impor¬ 
tant to use results from the same experiment in evalu¬ 
ating the right-hand side of Eq. ( [Tl] ). We use cumulant 
results from CMS [i]. CMS has not published U2{6}, 
but ATLAS has observed [ 7 ] that U2{6} ~ U2{4} for all 
centralities, therefore we assume U2{6} = U2{4}. 

The response coefficients in the left-hand side of 
Eq. (Ill are calculated using hydrodynamics. The calcu¬ 
lation shown in Fig. is the same as in Ref. |TH]. It uses 
as initial condition a symmetric Gaussian density prohle, 
where the normalization is adjusted to fit the measured 
multiplicity dNch/dy of Pb-Pb collisions at the LHC in 
the corresponding centrality class. This symmetric pro¬ 
file is deformed in order to produce anisotropic flow in 
the desired harmonic}^ We assume uniform longitudinal 
expansion [ 32 ] • With these initial conditions, we solve 


® For instance, X4 obtained by introducing an elliptic deforma¬ 
tion and calculating X 4 = 774 /( 1 , 2 )^. 


ideal hydrodynamics or second order viscous hydrody¬ 
namics | 43 | with constant shear viscosity over entropy 
ratio rj/s = 0.08 [ 33 ]. The equation of state is taken 
from Lattice QCD | 35 |. The initial time of the calcula¬ 
tion is To = 1 fm/c and the freeze-out temperature | 46 j 
is Tfo = 150 MeV. Anisotropic flow, u„, is calculated 
at freeze-out. It is averaged over particles in the inter¬ 
val Pt > 0.3 GeV/c, corresponding to the CMS accep¬ 
tance |3]. 

Figure[^shows that hydrodynamics naturally captures 
the sign, the magnitude, and the centrality dependence 
of all four nonlinear response coefficients. Experimental 
results differ from hydrodynamic calculations only for the 
most central bins | 41 ] , where the linear part typically be¬ 
comes larger than the nonlinear part and their correlation 
can no longer be neglected. 

The order of magnitude of the hydrodynamic result 
can be understood simply. At fixed, large px, ideal hy¬ 
drodynamics predicts HSIIIH] X 4 = 5, X 5 = 1 , X 62 = g, 
X63 = 5, X7 = 5- However, after averaging over px, X4 is 
multiplied by {V2)/{v2)‘^ > 1, where brackets now denote 
an average over px in a single hydro event. This is the 
reason why the results shown in Fig. are larger than 
the fixed-pT prediction. Since V2 and V3 have similar 
Px dependences, the enhancement factor is roughly the 
same for all quadratic response terms: panels (a), (b), 
(d) show that xs ~ 2^4, xes X4, in agreement with 
the above values. The enhancement from averaging over 
Px is larger for cubic response terms than for quadratic 
terms, but it is similar for both cubic terms: panels (c) 
and (e) show that X7 3 X 62 , also in agreement with the 
above values. 

A full hydrodynamical calculation gives results which 
differ somewhat from the naive predictions above. Coef¬ 
ficients from ideal hydrodynamics have a slight centrality 
dependence which is not captured by these formulas m- 
Viscous hydrodynamics predicts lower coefficients than 
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FIG. 3. (Color online) n5{4>23}, V 6 {'^ 2 }, r’el'J'a} and n7{4>23}, 
averaged over charged particles with pr > 0.3 GeV/c, as a 
function of centrality in Pb-Pb collisions at 2.76 TeV. vela's} 
has been shifted to the right by 1 % for sake of clarity. 


ideal hydrodynamics. Effect of viscosity, however, cancel 
to a large extent in the ratios: they are much smaller on 
X 4 = va/{v 2 Y than on and ( 1 ^ 2 )^ individually [T5] . 

Nonlinear response coefficients are mostly determined 
at freeze-out [18] , which is probably the least understood 
part of hydrodynamic calculations. While they depend 
little on the details of the initial profile or of the hydro- 
dynamic evolution, they depend rather strongly on the 
freeze-out temperature [47] . Similarly, the dependence of 
our results on viscosity is mostly through the viscous cor¬ 
rection to the momentum distribution at freeze-out |48| . 
The momentum distribution at freeze-out is not con¬ 
strained theoretically [49l EOj , and the quadratic ansatz 
used in this calculation is not favored by previous studies 
of U 4 m- Our calculation does not involve bulk viscos¬ 
ity, which is likely to be important at freeze-out [52fl55j . 
Finally, our results are quite sensitive to the value of the 
freeze-out temperature. Further studies are needed in or¬ 
der to pin down the sensitivity of response coefficients to 
model parameters. 


IV. PREDICTIONS 


Figure displays U5{4'23}, U 6 { 4 ' 2 }, and 

w{^ 23 }- Out of these four quantities, only U 6 { 4 ' 2 } has 
been measured by CMS |1|. We predict U 5 {'I' 23 } and 
using Eqs. (|^ and (|^, where we take U 5 { 4 > 5 } 
and ueldle} from CMS |4] and P 235 and psQ from AT¬ 
LAS |5|. 

Finally, W 7 {'I' 23 } is obtained from the last line of 
Eq. We use the viscous hydrodynamic calcula¬ 

tion for X 7 shown in Fig. [2 (e). We again assume that 
the magnitudes of V 2 and V 3 are independent, that is, 
(|k 2 |^|V 3 ^|) ~ ( 1141 ^) ( 103 ^ 1 ), and we estimate the moments 


using Eq. (12) and CMS data |3]. We anticipate that the 
absolute experimental error on 1:7 { 4123 } should be simi¬ 
lar to the error on U6{4'2}. This error is of the order of 
0.01%. The predicted values of U 7 {'I' 23 } is 0.05% in the 
25-30% centrality range, larger than the error. We there¬ 
fore expect that a nontrivial W 7 {'I' 23 } could be measured 
in midcentral Pb-Pb collisions at the LHC. 


V. CONCLUSION 

Harmonics U 4 and higher can be measured either with 
respect to their own planes or with respect to lower har¬ 
monic planes. We have clarified the relation between 
these projected harmonics and the so-called event-plane 
correlations. 

We have shown that under fairly general assumptions, 
measurements of higher harmonics can be combined with 
measurements of V 2 and V 3 in a way that eliminates the 
dependence on the initial state, and can be directly com¬ 
pared with hydrodynamic calculations. Experimental re¬ 
sults for V 4 , U 5 and Vq are in good agreement with vis¬ 
cous hydrodynamic calculations. We have argued that 
vr could be measured, and presented quantitative pre¬ 
dictions. 

On the experimental side, analyses should be repeated 
using the scalar-product method, whose result may differ 
significantly from the event-plane method for higher har¬ 
monics m- On the theoretical side, we hope that stud¬ 
ies of higher harmonics will help constrain the theoretical 
description of the fluid close the freeze-out temperature, 
which is poorly understood at present. 
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Appendix A: Analysis 

The flow observables in Eq. §, (§ and ([^ are ex¬ 
pressed in terms of moments of the distribution of I 4 ,. A 
generic moment is of the form [31j 



where and are integers and azimuthal symmetry 
implies Y^„nkn = For instance, (V 4 (V' 2 *)^) 

corresponds to = 1 , I 2 = 2 ; (IE 2 I®) corresponds to 
k 2 = h = 3; corresponds to k 2 = h = 2, 

ks = I3 = 1. 

We now describe a simple procedure for measuring 
these moments |31j , which generalizes the scalar-product 
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method [SSj. We define in each collision the flow vector 
by 


(A2) 

3 


where the sum runs over N particles seen in a refer¬ 
ence detector, and (pj are their azimuthal angles. One 
measures Qn in two different parts of the detector 
(“subevents”) A and which are symmetric around 
midrapidity and separated by a gap in pseudorapidity in 
order to suppress nonflow correlations [T31[S71[^. The 
moment (A11 is then given by 


Applied to Eq. ([^, this gives: 

^ I = ^(^{Q5aQ2bQ*3b) 


(A4) 


3B/ 


The scalar-product method thus uses the magnitude of 
the flow vector while the traditional event-plane 
method [59| only uses its azimuthal angle. One can 
symmetrize the numerator of Eq. (A4) over A and B 


to decrease the statistical error. Instead of 2 symmetric 
subevents, one can use 3 non-symmetric subevents, as 
described in Ref. m- 

Finally, analyses can be done differentially (in px bins, 
for identified part icles , etc.). For the differential analysis, 
one replaces Eq. (A4) by: 


t’5{'I'23} = 


Re{e^^‘^Q*BQ*3B) 

\/Rb{Q2aQ3aQ2bQ3b) 


(A5) 


M = 



{QnA)'‘-iQ*r.B) 



where the average in the numerator is now an average 
over particles in the considered bin, with azimuthal angle 
(A3) ip, instead of an average over events. 
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